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Abstract 

The Lane-Emden equation has been used to model several phenomenas in theoretical physics, math¬ 
ematical physics and astrophysics such as the theory of stellar structure. This study is an attempt to 
utilize the collocation method with the Rational Chebyshev of Second Kind function (RSC) to solve 
the Lane-Emden equation over the semi-infinit interval [0,+oo). According to well-known results and 
comparing with previous methods, it can be said that this method is efficient and applicable. 
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1 introduction 


Solving science and engineering problems appeared in unbounded domains, scientists have studied about it 
in the last decade. Spectral method is one of the famous solution. Some researchers have proposed several 
spectral method to solve these kinds of problems. Using spectral method related to orthogonal systems in 
unbounded domains such as the Hermite spectral method and the Laguerre Method is called direct approx¬ 
imation m- 

In an indirect approximation, the original problem mapping in an unbounded domain is changed to a problem 
in a bounded domain. Gua has used this method by choosing a suitable Jacobi polynomials to approximate 
the results Gann- There is Another method called domain truncation. In this method, the [0,+oo) to [0,L] 
and (—oo,+oo) to [-L,+L] are replaced by choosing a sufficient large L [l^j. Using spectral method based on 
the rational approximation is one of the direct approximation has been used in this paper. Authors of [IMS] 
have developed some spectral methods on unbounded intervals by using some mutually orthogonal systems of 
rational functions. Boyd m has proposed a new spectral basis named Rational Chebyshev functions on the 
semi-infinite interval by mapping to the Chebyshev polynomials. Author of [TS] has suggested and discussed 
a set of Legendre rational functions which are mutually orthogonal in L 2 (0, +oo) with a non-uniform weight 
function uj(x) = {x + l)~ 2 . Authors of JT7H2T] have applied some Rational Chebyshev and Legendre functions 
by tau and collocation method which were used to solve nonlinear ordinary differential equations on semi¬ 
infinite intervals. Author of mm have obtained the operational matrices of the derivative and product of 
rational Chebyshev and Legendre functions to reduce the solution of these problems to the solution of some 
algebraic equations. 


The Lane-Emden equation has been used to model several phenomenas in theoretical physics, mathe¬ 
matical physics and astrophysics as the theory of stellar structure. It appears in other contexts, for example 
in case of radiatively cooling, self-gravitating gas clouds, in the mean-field treatment of a phase transition 
in critical absorption or in the modeling of clusters of galaxies. It has been classified as singular initial 
value problems in the second order ordinary differential equations. This equation is to honor astrophysicists 
Jonathan Homer Lane and Robert Emden which is called Lane-Emden. [221127] 

Consider the Poisson equation and the condition for hydrostatic equilibrium [28l ~ l30l (42] : 


dP 

dr 


= ~P 


GM(r ) 
^2 ’ 


(1) 


d,M (7') 
dr 


= 4?rpr 2 


( 2 ) 


where G is the gravitational constant, P is the pressure, M(r) is the mass of a star at a certain radius r, and p 
is the density, at a distance r from the center of a spherical star [21].The combination of these equations yields 
the following equation, which as should be noted, is an equivalent form of the Poisson equation j28ll30[ (42j: 


1 d f r 2 dP\ 

I s dr (.717J = ~ 4 ’ rGp ' 


(3) 


From these equations one can obtain the Lane-Emden type equation through the simple assumption that 
the density is simply related to the density, while remaining independent of the temperature. We already 
know that in the case of a degenerate electron gas, the pressure and density are p ~ Ps, assuming that such 
a relation exists for other states of the star, we are led to consider a relation of the following form [28] : 


P = Kp 1+ ™, (4) 

where K and m are constants and m is the polytropic index which is related to the ratio of specific 
heats of the gas comprising the star.Based upon these assumptions we can insert this relation into our first 
equation for the hydrostatic equilibrium condition and from this equation |28H301 142] we have 


K (to + 1) ' u__i 

4t tG 


1 d 

( r 2^y\ 

r 2 dr 

K dr) 


= -y 


(5) 


where the additional alteration to the expression for density has been inserted with A representing the central 
density of the star and y that of a related dimensionless quantity that are both related to p through the 
following relation P%j - fo0l 1421 

P = A 2/ m (6) 
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Additionally, if place this result into the Poisson equation, we obtain a differential equation for the mass, 
with a dependence upon the polytropic index to. Though the differential equation is seemingly difficult to 
solve, this problem can be partially alleviated by the introduction of an additional dimensionless variable x 
,given by the following: 

r = ax (7) 

( 8 ) 


a 


A (to + 1) i 


4t tG 


A- 


--1 


Inserting these relations into our previous equations we obtain the famous form of the Lane-Emden type 
equations, given in the following: 


x 2 d'X \ dx ) 


(9) 


Taking these simple relations we will have the standard Lane-Emden equation with g(y) = y m [28H301 I42j . 


y" + -y' + y m = 0, x > 0. 

X 


( 10 ) 


At this point it is also important to introduce the boundary conditions which are based upon the following 
boundary conditions for hydrostatic equilibrium and normalization consideration of the newly introduced 
quantities x and y. What follows for r = 0 is 


r = 0—>£ = 0,p = A—>y(0) = l. (11) 

As a result an additional condition must be introduced in order to maintain the condition of equation(??) 
simultaneously: 

y' = o (12) 

In other words, the boundary conditions are as follows 

V (0) = 1, y' (0) = 0 (13) 

The values of to which are physically interesting, lie in the interval [0,5]. The main difficulty in the analysis 
of this type of equation is the singularity behavior occurring at x = 0. Exact solutions for equation (??) are 
known only form = 0,1 and 5. For other values of to the standard Lane-Emden equation is to be integrated 
numerically. Thus we decided to present a new and efficient technique to solve it numerically. 

This paper is organized as follows: In Section 2 we survey several methods that have been used to solve 
Lane-Emden type equations. In Section 3, the properties of Chebyshev Function of the Second Kind and the 
way to construct the collocation technique for this type of equation are described. In Section 4 the proposed 
method is applied to some types of Lane-Emden equations, and a comparison is made with the existing 
analytic or exact solutions that were reported in other published works in the literature. Finally we give a 
brief conclusion in the last section. 


2 The Methods have been proposed to solve Lane-Emden type 
equations 

Many analytical methods have been applied to solve Lane-Emden type equations. The major problem has 
been raised in the singularity of the equations at x = 0. 

Bender et al.[22] have shown the perturbative technique by using d-method has given excellent result 
when applied to difficult non-linear differential equation such as Lane-Emden. 

Kara et al.[52] have shown the Lagrangian method to solve Emden equation. They have categorized in 
non-linearzable equations and one year later pT3j they have found general solutions of Lane-Emden equation 
by determining the Lie point symmetries and, hence, the Lie algebra it generates. 

Shawagfeh[M| has approached the problem differently by utilizing the newly developed Adomian de¬ 
composition method to solve Lane-Emden equation and has used this method to derive a nonperturbative 
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approximate solution in the form of power series which is easily computable and accurate. To accelerate the 
convergence he has used the Pade‘ approximates method. 

Mandelzweig et al. [35] have shown that the quasilinearization method (QLM) gives excellent results when 
applied to different nonlinear ordinary differential equations in physics, such as Lane-Emden equation. 

Wazwaz [3T has employed a reliable algorithm rests mainly on the Adomian decomposition method 
with an alternate framework designed to overcome the difficulty of the singular point and he has proposed 
framework has been applied to a generalization of Lane-Emden equations. 

Liao [35] has introduced an analytic algorithm logically contains the well-known Adomian decomposition 
method, provides for Lane-Emden equation. This algorithm itself provides us with a convenient way to adjust 
convergence regions even without Pade technique. 

He [361 has presented by the semi-inverse method, a variational rinciple is obtained for the Lane-Emden 
equation, which gives much numerical convenience when applying finite element methods or Ritz method. 

Ramos m has uttered Linearization methods provide piecewise linear ordinary differential equations 
which can be easily integrated, and provide accurate answers even for hypersingular potentials, for which 
perturbation methods diverge and he has shown the applicability and accuracy of linearization methods for 
initial-value problems in ordinary differential equations are verified on examples that include the Lane-Emden 
equation, and some other equations. 

Parand et al. |38| has expressed numerical method based on rational Legendre tau for solving the Lane- 
Emden nonlinear differential equation. They have provided the operational matrices of derivative and product 
of rational Legendre functions. 

Ramos [331 has employed Linearization methods result in linear constant-coefficients ordinary differential 
equations which can be integrated analytically, thus yielding piecewise analytical solutions and globally 
smooth solutions. 

Wazwaz m has shown the modified decomposition method which has been applied for analytic treatment 
of nonlinear differential equations such as Lane-Emden. This method accelerates the rapid convergence of the 
series solution, dramatically reduces the size of work, and provides the exact solution by using few iterations 
only without any need to the so-called Adomian polynomials. 

Yusefi [26] has presented a numerical method by using integral operator and convert Lane-Emden equa¬ 
tions to integral equations and then applying Legendre wavelet approximations. 

Ramos [41] has detailed series solutions of the Lane-Emden equation based on either a Volterra integral 
equation formulation or the expansion of the dependent variable in the original ordinary differential equation 
and has compared with series solutions obtained by means of integral or differential equations based on a 
transformation of the dependent variables. 

Aslanov |42| has constructed a recurrence relation for the components of the approximate solution and 
has investigated the convergence conditions of the Emden-Fowler type of equations. 

Razzaghi et al |43| have given numerical method based upon hybrid function approximations for solving 
nonlinear initial-value problems such as Lane-Emden type equations.The properties of hybrid of block-pulse 
functions and Lagrange interpolating polynomials have presented and and have utilized to reduce the com¬ 
putation of nonlinear initial-value problems to a system of non-algebraic equations. 

Dehghan et al.03] have investigated The Lane-Emden using the variational iteration method and it have 
been shown the efficiency and applicability of this procedure for solving this equation. 

Hashinr et al. [45] have obtained approximate and/or exact analytical solutions of the generalized Emden- 
Fowler type equations in the second-order ordinary differential equations (ODEs) by homotopy-perturbation 
method (HPM). 

Parand et al [JB] have represented a pseudospectral technique is to solve the Lane-Emden type equations 
on a semi-infinite domain. The method is based on rational Legendre functions and Gauss-Radau integration. 

Ramos [47] has shown that the appearance of noise terms in the decomposition method is related to both 
the differential equation and the manner in which the homotopy parameter is introduced, especially for the 
Lane-Emden equation. 

Bataineh et al [H] have stated approximate and/or exact analytical solutions of singular initial value 
problems (IVPs) of the Emden-Fowler type in the second-order ordinary differential equations (ODEs) are 
obtained by the homotopy analysis method (HAM). 

Ozi§ et al.|43] have declared, approximate-exact solutions of a class of Lane-Emden type singular IVPs 
problems, by the variational iteration method. 
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Singh et al. |50j have shown an efficient analytic algorithm for Lane-Emden type equations using modified 
homotopy analysis method, which is different from other analytic techniques as it itself provides us with a 
convenient way to adjust convergence regions even without Pade technique. 

Parand et al. m have provide a pseudo spectral method for Lane-Emden equation which based on some 
orthogonal functions. 

Parand et al. [53] have imparted a method based on the rational Legendre functions and Gauss-Radau 
integration to solve Lane-Emden type equations on a semi-infinite domain. 

Geng et al.|54] have investigated the nonlinear singular initial value problems including generalized Lane- 
Emden-type equations by combining homotopy perturbation method (HPM) and reproducing kernel Hilbert 
space method (RKHSM). 

Karimi et al.|55) have presented a numerical method which produces an approximate polynomial solution 
for solving Lane-Emden equations as singular initial value problems. They have used an integral operator to 
convert Lane-Emden equations into integral equations the convert the acquired integral equation into a power 
series and finally, transforming the power series into Pad series form, they have obtained an approximate 
polynomial of arbitrary order to solve Lane-Emden equations. 

Gorder |56j has applied the 5-expansion method to a transformed Lane-Emden equation then the results 
have been transformed back, and has been recovered analytical solutions to the Lane-Emden equation of the 
second kind in a special case. The rapid convergence of the method results in qualitatively accurate solutions 
in relatively few iterations, as we see when we compare the obtained analytical solutions to numerical results. 

Wazwaz et al. SB] have provided a comparison of the Adomian decomposition method (ADM) with the 
variational iteration method (VIM) for solving the Lane-Emden equations of the first and second kinds. 

Yzbasi [59] has mentioned a collocation method based on the Bessel polynomials is presented for the 
approximate solution of a class of the nonlinear Lane-Emden type equations. 

Bhrawy et al. [50] have conveyed a shifted Jacobi pseudo-spectral method JPSM, along with the Boubaker 
Polynomials Expansion Scheme (BPES), for solving a nonlinear Lane-Emden type equation. 

Bhrawy et al. m have related a shifted Jacobi-Gauss collocation spectral method for solving the nonlinear 
Lane-Emden type equation. 

Caglar et al. [52] have solved a time-dependent heat-like Lane-Emden equation by using a non-polynomial 
spline method. 

Shen [53] has combined the compactly supported radial basis function (RBF) collocation method and 
the scaling iterative algorithm to compute and visualize the multiple solutions of the Lane-Emden-Fowler 
equation on a bounded domain !!c R 2 with a homogeneous Dirichlet boundary condition. 

Pandey et al. (64| have developed an efficient numerical method for solving linear and nonlinear Lane- 
Emden type equations using Legendre operational matrix of differentiation. 

Pandey et al. [55] have applied the First Bernstein operational matrix of differentiation derived using 
Bernstein polynomials to solve the linear and nonlinear differential equations of Lane-Emden type. 

Biezuner et al. [67] have introduced an iterative method to compute the first eigenpair (A p , e p ) for the 
p-Laplacian operator with homogeneous Dirichlet data as the limit of (/i q ,u q ) as q —> p~, where u q is the 
positive solution of the sublinear Lane-Emden equation —A p u q = fiqU ^ -1 with the same boundary data. 

Boubaker et al. [53] have used the Boubaker Polynomials Expansion Scheme (BPES) in order to obtain 
analytical-numerical solutions to the Lane-Emden initial value problem of the first kind and second kind. 

Rismani et al. [59] have utilized the improved Legendre-spectral method to solve Lane-Emden type 
equations. 

Jalab et al. m have shown a numerical method based on neural network, for solving the Lane-Emden 
equations singular initial value problems. The numerical solution has been given for integer case and non 
integer case. The non integer case is taken in the sense of Riemann-Liouville operators. 

Wazwaz et al. eh have used the systematic Adomian decomposition method to handle the integral form 
of the Lane-Emden equations with initial values and boundary conditions and have confirmed their belief 
that the Adomian decomposition method (has provided) provides efficient algorithm for analytic approximate 
solutions of the equation. 

Karimi Vanani et al. m have presented a numerical algorithm based on an operational Tau method 
(OTM) for solving the Lane-Emden equations as singular initial value problems. 

Parand et al [73] have introduced the Bessel orthogonal functions as new basis for spectral methods and 
also present an efficient numerical algorithm based on them and collocation method to solve these well-known 
equations. 
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Kaur et al. m have provided a technique to investigate the solutions of generalized nonlinear singular 
Lane-Emden equations of first and second kinds by using a Haar wavelet quasi-linearization approach. 

Oztiirk et al. m have stated the numerical solution of Lane-Emden equations by using truncated shifted 
Chebyshev series together with the operational matrix. 

Razzaghi |78J has proposed a numerical method based on hybrid function approximations and Bernoulli 
polynomials has been presented and utilized to reduce the computation of nonlinear initial-value problems 
such as Lane-Emden type equations to a system of equations. 

Li et al CS3 have computed and visualized multiple solutions of the Lane-Emden equation on a square 
and a disc, using Legendre and Fourier-Legendre pseudospectral methods based on the Liapunov-Schmidt 
reduction and symmetry-breaking bifurcation theory. 

Wazwaz [8T] has established the Volterra integro-differential forms of the Lane-Emden equations and has 
used the variational iteration method (VIM) to effectively treat these forms. 

Rach et al. [82] have considered the coupled Lane-Emden boundary value problems in catalytic diffusion 
reactions by the Adomian decomposition method. 

Mohammadzadeh et al. has shown three numerical techniques based on cubic Hermite spline functions 
for the solution of Lane-Emden equation. 

Wazwaz et al. [84] have investigated systems of Volterra integral forms of the Lane-Emden equations 
have used the systematic Adomian decomposition method to handle these systems of integral forms. 

Oztiirk et al [55] have shown numerical method depends on collocation method and based on first taking 
the truncated Hermite series of the solution function, transforms Lane-Emden type equation and given 
conditions into a matrix equation,solving the system of algebraic equations using collocation points and 
found the coefficients of the truncated Hermite series. 

Wazwaz et al. [SB] have presented new alternate derivations for the Volterra integral forms of the Lane- 
Emden equation for the shape factor of k=l, where an alternate derivation for L’Hospital’s formula will be 
developed. 

Gurbiiz et al. m have mentioned the numerical method is based on the matrix relations of Laguerre 
polynomials and their derivatives, and reduces the solution of the Lane-Emden type functional differential 
equation to the solution of a matrix equation corresponding to system of algebraic equations with the unknown 
Laguerre coefficients. 


3 Properties of Chebyshev Function of the Second Kind 


The chebyshev polynomials in second kind are defined on [—1,1] by 

U„ (. x ) = sm ( n + 1 ^ , x = cos6>,n = 0,l,2,--- 
sin 0 


(14) 


with respect to weight function w ( x ) = y/\ — x 2 and can be determined following recurrence formula: 

Uq ( x ) = 1, U\ ( x ) = 2x 

U n (x) = 2xU n _ 1 (x) - U n —2 (x),n = 2,3,4, • • • (15) 

We can obtain another formula for derivative of U n (x) as 


d (n + 2)U n -i(x) — nU n+ i(x) 

dx U " (x) =- 2 ( 1 ^)- 


(16) 


Since the range [0, + 00 ) quite often more convenient to use than the range [—1,-f-l], we sometimes map 
the independent variable x in [0, + 00 ) to the variable s in [—1, +1] by the transformation 


x — 1 s + 1 

s = -- or x = 


x +1 


1 — s 


(17) 


and this leads to a shifted Chebyshev polynomial (of the second kind) U* of degree n in x on [0, + 00 ) given 
by 

U*(x) = U n (s) = U n (^\ (18) 
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with respect to weight function w*(x) = H§js- The U* has been called Rational Chebyshev of Second kind 
(RCS). From (??) and (??), we may deduce the recurrence relation for U* in the form 

' x — 1' 


K (x) = 2 


x +1 


with initial conditions 


Uq (x) = 1, Ut (x) = 2 

As same like (??) we can obtain derivative of U*(x) as follow 


mJ*(x) = „ 

dx (l + x) 


K-i (x) - U* n _ 2 (x) ,n = 2,3,4, 

x — 1 

x + 1 / 

2 (n + 2)U* 1 (x) — nU* +1 (x) 


2(1-(i^) ) 


(19) 

( 20 ) 

( 21 ) 


Theorem 1. In the Gauss Chebyshev formula 

\ n 

/ f{x)w(x)dx ~ Y, A k f(x k ) (22) 

_1 fe=l 

where{xk} are the n zeros of <f n (x),the coefficients A k are as follows: 
for w(x) = (1 — x)“ 2 ; cj) n (x) = U n (x ) and A k = (1 — x%) 

Proof. The complete proof has been shown by Mason et al. in [55] □ 

Let w*(x) = ^+ 3 ) a denote a non-negative, integrable, real-valued function over the interval I = [0, +oo). 
We have defined 

LIj, (/) = {«:/—>■ K|n is measurable and IHI^, < oo} , (23) 

where 




\v (x)|“w;* (x) dx 


(24) 


is the norm induced by the scalar product 

OO 

(u,v} w „ = J u(x)v(x)w* (x)dx 
o 

Thus {U* (x)} n>0 denote a system which are mutually orthogonal under (??), i.e 

(U*{x),U* m {x)) w , = ^5 nrn 

where S nm is the Kronecker delta function. This system is complete in l? w , (/). For any function u £ L (/) 
the following expansion holds 


(25) 


(26) 


+N 

/(x) = a i U i( X ) 

i=-N 

<«. Uf) w , 

a k = 


(27) 


(28) 


\m\c 

The a k s are the expansion coefficients associated with the family {17* (x)}. Now we can define an orthogonal 
projection based on the transformed Chebyshev functions in second kind as given below: Let 

U* N = span{Uo(x),Uf(x),...,U*(x)} 

The L 2 (J)-orthogonal projection : L 2 (I) —> is a mapping in a way that for any y £ L 2 (/), 

(iNy-y,I^ = O,V0e 13* N 

or equivalently, 

N 

£ N y = J2 a i u *(x) (29) 


i=0 
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4 Applications 

In this section we have applied Rational Chebyshev of Second kind Collocation method for the computation 
of Lane-Emden type equations based on the Rational Chebyshev of Second kind functions. In general the 
Lane-Emden type equations are formulated as 

ry 

y” + -y + f(x)g(x) = h(x), ax>0. (30) 


with initial conditions 

2/(0) = A, y'( 0) = B 

where a, A and B are real constants and fix), g(y) and h{x) are some given functions. 

We apply the Rational Chebyshev of Second kind functions collocation method to solve some well-known 
Lane-Emden type equations for various f(x), g(y), A and B, in two cases homogeneous (h{x) = 0) and 
non-homogeneous (h{x) ^ 0). 

To satisfy boundary conditions we can add the (??) by P(x) = A + Bx, therefore the (??) has been 
changed as follows: 


£nv{x) = P{x) + x 2 £ N y{x). 

Now for boundary conditions we have £nv(x) = A and yzf,Ny{x) = B when x tends to zero. 

To apply the collocation method, we have constructed the residual function by substituting y{x) by 
^N'y(x) in the Lane-Emden type equation (??): 

Res(:r) = J^6v2/(a0 + ~^A N y(x) + f{x)g{£ N y{x)) - h(x). 


There is no limitation to choose the point in collocation method. In this case, some points with equal dis¬ 
tance from each other in [0, q\ have been chosen. The amount of q has been provided by physical view of 
Lane-Emden equation. N + 1 points with equal distance from each other in [0, q] have been substituted in 
Res{x ): 


Res(xj) = 0, j = 0,1,2,3,..., N 

A non-linear system of equations has been obtained. Some numerical methods should be used as Newton’s 
method. The Maple version 17 has been used to solve the non-linear system of equations. Solving these 
types of equations system, Maple software has been used the Newton’s method and advanced algorithm. 
The Maple /solve has been used for all of unknown coefficients with an initial zero value. The unknown 
coefficients cq will be obtained and approximation of ^N'yix) is the result. 


4.1 Example 1 (The standard Lane-Emden equation) 

According to Eq.(??), if f(x) = l,g(y) = y m , A = 1 and B = 0 the standard Lane-Emden equation has been 
defined as follows: 

y"{x) + ^y'(x) + y m (x) = 0,x ^ 0 (31) 

with boundary conditions 


2 /( 0 ) = 1 , 2 /( 0 ) = 0 


where m > 0 is constant. 

In Eq. (??) for m = 0 exact solution is y{x) = 1 — ^ x 2 , for m = 1 is y(x) = 
y(x) = (l + 2 - 

The Rational Chebyshev of Second kind Collocation method has been applied to 


Sln( - a, ' ) and for m = 5 is 

X 


solve Eq. (??) for m = 



1.5, 2,2.5, 3 and 4 because there is no analytic exact solution for other case of to values so the residual function 
has been constructed as: 


d? ~ 2d- /- \ m 

Res(x) = ^ ^nv(x) + ~^nv{x) + \^Ny{x)J . (32) 

As it was mentioned, N+l points with equal distance from each other in [0, q] have been substituted in Res(x): 

Res(xj) = 0, j = 0,1,2,3,..., N 

A non-linear system of equations has been obtained. Some numerical methods should be used as Newton’s 
method. The Maple version 17 has been used to solve the non-linear system of equations. Solving these 
types of equations system, Maple software has been used the Newton’s method and advanced algorithm. 
The Maple f solve has been used for all of unknown coefficients with an initial zero value. The unknown 
coefficients cii will be obtained and approximation of ry(x) is the result. 

Table [T] have shown comparison of the first zeros of standard Lane-Emden equations, for the current 
method and exact numerical values given by Horedt [85] for m = 1.5, 2, 2.5,3 and 4 with N = 20 ???? have 
represented approximations of y(x) for the standard Lane-Emden equation for m = 2, 3 respectively obtained 
by Horedt [82] and the method proposed in this paper 

and it shows the accuracy from this method and the exact value is equal up to 8 digits. 

?????? have represented approximations of y{x) for the standard Lane-Emden equation for m = 1.5, 2.5, 3.5 
and 4 respectively obtained by by Horedt [55] and the method proposed in this paper 

?? has shown the coefficients of the Chebyshev of the Second kinds functions obtained by the curent 
method for to = 1.5,2, 2.5, 3 and 4 of the standard Lane-Emden equation. 

The result graph of the standard Lane-Emden equation for to = 1.5, 2,2.5, 3, 3.5 and 4 has shown in ??. The 
logarithmic graph of absolute coefficients of Chebyshev of the Second kinds of functions of standard Lane- 
Emden equation for m = 1.5 is shown in ??. The coefficients in ?? has shown that there is an appropriate 
convergence rate for the new method. 

Table 1: Comparison of the first zeros of standard Lane-Emden equations, for the current method and exact 
numerical values given by Horedt |89j 


m 

N 

Present method 

Exact Value 

1.50 

20 

3.65375374E + 00 

3.65375374E + 00 

2.00 

20 

4.35287460E + 00 

4.35287460E + 00 

2.50 

20 

5.35527546E + 00 

5.35527546E + 00 

3.00 

20 

6.89684862E + 00 

6.89684862E + 00 

4.00 

20 

1.49715463E + 01 

1.49715463E + 01 


Table 2: Comparison of y{x) values of standard Lane-Emden equation, for the current method and exact 
values given by Horedt for m=1.5 [55] 


X 

Present method 

Exact Value 

Error 

0.00 

0.10 

0.50 

1.00 

3.00 

3.60 

3.65 

3.6537537 

1.00000000E + 00 
9.98334583E - 01 
9.59103857E - 01 
8.45169755E - 01 
1.58850614E - 01 
1.10779151E - 02 

7.62752205E - 04 
2.44970506E - 10 

1.00000000E + 00 
9.983346007? - 01 
9.591039007? - 01 
8.451698007? - 01 
1.588576007? - 01 
1.109099007? - 02 
7.639242007? - 04 
0.000000007? + 00 

0.000000007? + 00 
1.709599807? - 08 
4.306351907? - 08 
4.451393207? - 08 
6.986144987? - 06 
1.307488477? - 05 
1.171994717? - 06 
0.000000007? + 00 
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Table 3: Comparison of y{x) values of standard Lane-Emden equation, for the current method and exact 
values given by Horedt for m=2 [89] 


X 

Present method 

Exact Value 

Error 

0.00 

0.10 

0.50 

1.00 

3.00 

4.00 

4.30 

4.35 

4.35287460 

1.00000000E + 00 
9.98334988 E - 01 
9.593527057? - 01 
8.486541037? - 01 
2.418240787? - 01 
4.884013237? - 02 

6.810928617? - 03 
3.660296767? - 04 

o.ooooooooe; + oo 

1.00000000.E + 00 
9.983350007? - 01 
9.593527007? - 01 
8.486541007? - 01 
2.418241007? - 01 

4.884015007? - 02 
6.810943007? - 03 
3.660302007? - 04 

o.ooooooooe; + oo 

o.ooooooooe; + oo 

1.175889507? - 08 
4.980445067? - 09 
2.862989077? - 09 
2.169664407? - 08 
1.774841127? - 08 

1.438785257? - 08 
5.243432657? - 10 
0.000000007? + 00 


Table 4: Comparison of y{x ) values of standard Lane-Emden equation, for the current method and exact 
values given by Horedt for m=2.5 [89] 


X 

Present method 

Exact Value 

Error 

0.00 

0.10 

0.50 

1.00 

4.00 

5.00 

5.30 

5.355 

5.35527546 

1.00000000.E + 00 
9.982001617? — 01 
9.59597754 E - 01 
8.519441997? - 01 

1.376807517? — 01 
2.901986397? - 02 
4.259864197? - 03 
2.101105377? — 05 

o.ooooooooe; + oo 

1.00000000E + 00 
9.983354007? - 01 
9.595978007? - 01 
8.519442007? - 01 

1.376807007? - 01 
2.901919007? - 02 
4.259544007? - 03 
2.100890007? — 05 
0.000000007? + 00 

o.ooooooooe; + oo 

1.352389517? - 04 
4.553368097? - 08 
8.719049937? - 10 
5.093510707? - 08 
6.738935827? - 07 
3.201856447? - 07 
2.153671487? - 09 

o.ooooooooe; + oo 


Table 5: Comparison of y{x) values of standard Lane-Emden equation, for the current method and exact 
values given by Horedt for m=3 [89] 


X 

Present method 

Exact Value 

Error 

0.00 

0.10 

0.50 

1.00 

5.00 

6.00 

6.80 

6.90 

1.00000000# +00 
9.983358207? - 01 
9.598390607? - 01 
8.550575607? - 01 
1.108198307? — 01 
4.373798007? - 02 
4.167780007? — 03 
— 1.336500007? — 04 

1.000000007? + 00 
9.983358007? - 01 
9.598391007? - 01 
8.550576007? - 01 
1.108198007? - 01 
4.373800007? - 02 
4.167800007? - 03 
3.600000007? - 05 

0.000000007? + 00 
2.00000000# - 08 
4.000000007? - 08 
4.00000000# - 08 
3.000000007? - 08 
2.000000007? - 08 
2.00000000# - 08 
1.696500007? — 04 


Table 6: Comparison of y{x) values of standard Lane-Emden equation, for the current method and exact 
values given by Horedt for m=4 [39] 


X 

Present method 

Exact Value 

Error 

0.00 

0.10 

0.20 

0.50 

1.00 

5.00 

10.00 

14.00 

14.90 

1 . 00000000 # + oo 

1.015850887? + 00 
9.924250327? - 01 

9.603108917? — 01 
8.608138127? - 01 
2.359227317? - 01 

5.967377177? - 02 
8.332930357? - 03 
5.766738367? - 04 

1.00000000#+ 00 
9.983367007? - 01 
9.933862007? - 01 
9.603109007? - 01 
8.608138007? - 01 
2.359227007? - 01 
5.967274007? - 02 
8.330527007? - 03 
5.764189007? - 04 

o.oooooooo# + oo 

1.751418157? — 02 

9.611679057? - 04 
8.819220027? - 09 
1.220797517? — 08 
3.106852907? — 08 
1.031723177? - 06 
2.403346997? - 06 
2.549359167? - 07 
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Table 7: Coefficients of the RCS functions of the standard Lane-Emden equations for to = 1.5, 2, 3 and 4 


i 

Cti 


m = 1.5 

m = 2 

m = 2.5 

in = 3 

m = 4 

0 

-8.35020235E-02 

-7.37260675E-02 

-5.84134154E-02 

-6.00541969E-01 

-1.80533469E-02 

1 

4.14892772E-02 

4.61551717E-02 

4.24877591E-02 

1.07145455E+00 

1.60801870E-03 

2 

1.06527872E-02 

-5.59588313E-04 

9.68579002E-03 

-1.38004834E+00 

4.91151095E-02 

3 

-6.46059281E-03 

-2.88264289E-04 

-5.69688651E-03 

1.57447823E+00 

-6.29427069E-02 

4 

-9.23628705E-03 

-2.03783735E-02 

-2.83354601E-03 

-1.61608168E+00 

5.51425972E-02 

5 

-5.14378683E-03 

2.77045133E-03 

-1.88446376E-03 

1.49866216E+00 

-5.60724094E-02 

6 

-3.99197152E-04 

-1.30872723E-02 

2.47503817E-03 

-1.28488320E+00 

5.00091630E-02 

7 

2.35669147E-03 

6.02895464E-03 

4.20009874E-04 

1.01947325E+00 

-4.03638627E-02 

8 

3.03971126E-03 

-8.47557449E-03 

1.86252286E-03 

-7.46940934E-01 

3.09682458E-02 

9 

2.53242065E-03 

4.00096435E-03 

-1.11760997E-04 

5.05881970E-01 

-2.19137068E-02 

10 

1.67607539E-03 

-4.94521907E-03 

7.92650268E-04 

-3.15031208E-01 

1.42286699E-02 

11 

9.36473928E-04 

1.94713984E-03 

-1.16018023E-04 

1.79413013E-01 

-8.51801989E-03 

12 

4.51996582E-04 

-2.06474187E-03 

3.37002942E-04 

-9.27123693E-02 

4.63985922E-03 

13 

1.89757483E-04 

7.66054296E-04 

-1.36758208E-05 

4.30569575E-02 

-2.27214677E-03 

14 

6.90455024E-05 

-5.75592272E-04 

1.09800992E-04 

-1.76863574E-02 

9.90326051E-04 

15 

2.15040362E-05 

2.13683933E-04 

4.50594242E-06 

6.32155250E-03 

-3.76662896E-04 

16 

5.58921881E-06 

-9.82199992E-05 

2.08231737E-05 

-1.90188385E-03 

1.20837515E-04 

17 

1.16162729E-06 

3.53103253E-05 

1.59326696E-06 

4.63035504E-04 

-3.14970993E-05 

18 

1.75835253E-07 

-7.82964769E-06 

1.71120723E-06 

-8.26242326E-05 

6.09225231E-06 

19 

1.59423494E-08 

2.50723737E-06 

1.81759246E-07 

9.25997591E-06 

-7.06745003E-07 


Figure 1: Logarithmic graph of absolute coefficients |a.^| of RCS function of standard Lane-Emden equation 
form=1.5 
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Figure 2: Graph of standard Lane-Emden equation for m = 1.5, 2, 2.5,3 and 4. 


4.2 Example 2 (The isothermal gas spheres equation) 

According to Eq.(??), if f{x) = l,g(y) = e v , A = 0 and B = 0 the isothermal gas spheres equation has been 
defined as follows: 


2 

y"(x) H— y'(x) + e y (x) = 0, x ^ 0 
x 

with boundary conditions 


(33) 


2 /( 0 ) = 0, y'(0) = 0 

Davis (27) has discussed about ?? that can be used to view the isothermal gas spheres, where the temperature 
remains constant. 

A series solution have investigated by Wazwaz [31],Liao [35],Singh et al. [5D3 and Ramos m by using ADM, 
ADM, MHAM and series expansion, respectively: 


, V 1 2 ,1 4 8 

+ Mi 1 “ 2L6f 


122 
81.8! : 


61.67 

459.10!* 


„10 


(34) 


We have applied Chebyshev of Second kind Collocation method to solve Eq. (??) therefore we have con¬ 
structed the residual function: 


Res (a;) 


dx 2 


^NVix) 


2 d 
x dx 


InV(x) 


J,nv(x) 


(35) 


As it was mentioned, IV+1 points with equal distance from each other in [0, q] have been substituted in Res(x): 


Res(xj) = 0, j = 0,1,2,3,..., N 

A non-linear system of equations has been obtained. Some numerical methods should be used as Newton’s 
method. The Maple version 17 has been used to solve the non-linear system of equations. Solving these 
types of equations system, Maple software has been used the Newton’s method and advanced algorithm. 
The Maple /solve has been used for all of unknown coefficients with an initial zero value. The unknown 
coefficients cq will be obtained and approximation of £,nv(x) is the result. 
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?? has shown the comparison of y{x ) obtained by the method proposed in this paper with N = 40 and 
those obtained by Wazwaz [33]. As ADM is based on Taylor series ,it can be contributed that there is a 
very high level of accuracy around zero point. This accuracy will decrease while the distance of points starts 
to raise from zero point. The resulting graph of the isothermal gas spheres equation in comparison to the 
presented method and those obtained by Wazwaz[331 has been shown in ??. The logarithmic graph of the 
absolute coefficients of Chebysliev of Second kind functions of the standard isothermal gas spheres is shown 
in ??. This graph shows that the new method has a proper convergence rate. 

Table 8: Comparison of y{x), between current method and series solution given by Wazwaz [34] for isothermal 
gas sphere equation. 


X 

Present method 

Wazwaz 

Error 

0.00 

0.10 

0.20 

0.50 

1.00 

1.50 
2.00 

2.50 

0.00000000E + 00 
-1.66583386E-03 
-6.65336710E - 03 
—4.115395737? — 02 
— 1.588276787? — 01 
—3.380194257? - 01 
—5.598230047? - 01 
—8.063408717? - 01 

o.ooooooooe + oo 

— 1.665833907? - 03 
—6.653367107? - 03 
—4.115395687? - 02 

— 1.588273547? - 01 
—3.380131107? — 01 
—5.599626607? - 01 
—8.100196717? — 01 

0.000000007? + 00 
3.793941987? - 11 
4.231398227? - 13 

4.927138007? — 10 
3.238243947? - 07 
6.314460807? - 06 
1.396557647? - 04 
3.678800717? - 03 


Figure 3: Graph of isothermal gas sphere equation in comparison with Wazwaz solution [34] 
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Figure 4: Logarithmic graph of absolute coefficients \ai\ of RCS function of isothermal gas sphere equation. 


4.3 Example 3 

According to Eq.(??), if f{x) = 1, g{y) = sinh(y), A = 1 and B = 0 the equation has been defined as follows: 


2 

y"(x) + —y'{x) + sinh(y) = 0, x ^ 0 


(36) 


with boundary conditions 


1/(0) = W(0) = 0 

A series solution have investigated by Wazwaz [31] by using Adomian Decomposition Method (ADM) is: 

, . „ (e 2 — l)x 2 1 (e 4 — l)a; 4 1 (2e 6 + 3e 2 - 3e 4 - 2)x 6 1 (61e 8 + 104e 6 - 104e 2 - 61):r 8 

V {X)~ l~ V +— -- 7T 1 -77XXTX --5- ~ + ~ r ^ 


12e 480 e 2 


30240 


26127360 


we have purposed to RCS method to solve ?? therefore, we construct the residual function as follows: 

Res(z) = ^ NV{X) + ~c^ Ny ^ + sinh ^ Ny W>- 


(37) 


(38) 


N + 1 points with equal distance from each other in [0, q] have been substituted in R.es(x): 

Res(xj) = 0, j = 0,1,2,3 ,N 

A non-linear system of equations has been obtained. Some numerical methods should be used as Newton’s 
method. The Maple version 17 has been used to solve the non-linear system of equations. Solving these 
types of equations system, Maple software has been used the Newton’s method and advanced algorithm. 
The Maple /solve has been used for all of unknown coefficients with an initial zero value. The unknown 
coefficients will be obtained and approximation of £,Ny{x) is the result. 

?? has shown the comparison of y(x) obtained by the new method proposed in this paper with N=20, 
and those obtained by Wazwaz ]3!].The resulting graph of ?? in comparison to the presented method and 
those obtained by Wazwaz [35 are shown in ??. 
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Table 9: Comparison of y(x), between current method and series solution given by Wazwaz j.'jTl for Example 

3. 


X 

Present method 

Wazwaz 

Error 

0.00 

0.10 

0.20 

0.50 

1.00 

1.50 

2.00 

0.00000000E + 00 
9.98042842E - 01 

9.92189436E - 01 
9.51961094E - 01 
8.18242929E - 01 

6.25438765E - 01 
4.06623301E - 01 

0.00000000E + 00 
9.98042841E - 01 

9.92189435E - 01 
9.51961102E - 01 
8.18251667E - 01 
6.25891608E - 01 
4.13669104E - 01 

0.00000000E + 00 
9.80743042E - 10 
1.02811104E-09 

8.14090706E - 09 
8.73755808E - 06 
4.528431557? - 04 

7.045803127? - 03 


Figure 5: Graph of equation of Example 3 in comparing the presented method and Wazwaz solution 


4.4 Example 4 

According to Eq.(??), if f(x) = 1 ,g(y) = sin(y), A = 1 and B = 0 the equation has been defined as follows: 


2 

y”{x) + -y'{x) + sin(y) = 0, x > 0 


(39) 


with boundary conditions 


2 /( 0 ) = 1 , 2 /( 0 ) = 0 

A series solution have investigated by Wazwaz [M] by using Adomian Decomposition Method (ADM) is: 

1 


y(x ) ~ 1- kix 2 H- kikox 4 + k\ ( -fc?- kn ) x 6 + k\kn ( - k? - 

’ 6 120 1 \ 3024 1 5040 2 V 3265920 1 362880 2 


kn I a: 8 


+ k\ 


1781 


-kiki - 


1 


-kn ~ 


19 


82892800 1 2 399168000 2 23950080 


4 „.io 


kf ) x 


(40) 


where ki = sin(l) and &2 = cos(l) we have purposed to RCS method to solve ?? therefore, we construct the 
residual function as follows: 


Res(a:) = -^-y|jvy(a;) + --^-fiv2/(®) + sin (£ N y(x)). 
ax z x ax 


(41) 


As it was mentioned, N + 1 points with equal distance from each other in [0, q] have been substituted in 
Res(x ): 


Res(xj ) = 0, j = 0,1,2,3,..., N 
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A non-linear system of equations has been obtained. Some numerical methods should be used as Newton’s 
method. The Maple version 17 has been used to solve the non-linear system of equations. Solving these 
types of equations system, Maple software has been used the Newton’s method and advanced algorithm. 
The Maple /solve has been used for all of unknown coefficients with an initial zero value. The unknown 
coefficients cq will be obtained and approximation of £a ry{x) is the result. 

?? has shown the comparison of y(x ) obtained by the new method proposed in this paper with N=20, 
and those obtained by Wazwaz |M|.The resulting graph of ?? in comparison to the current method and those 
obtained by Wazwaz [34] are shown in ??. 

Table 10: Comparison of y(x), between current method and series solution given by Wazwaz [341 for Example 

4. 


X 

Present method 

Wazwaz 

Error 

0.00 

0.10 

0.20 

0.50 

1.00 

1.50 

2.00 

1.00000000#+ 00 
9.98597930# - 01 
9.94396268# - 01 
9.65177784# - 01 
8.63681129# - 01 
7.05045237# - 01 
5.06464502# - 01 

1.00000000#+ 00 
9.98597936# - 01 
9.94396273# - 01 
9.65177789# - 01 
8.63681103#-01 
7.05041925# - 01 
5.06372033# - 01 

0.00000000# + 00 
5.40042000# - 09 
5.12808607# - 09 
5.03923003# - 09 
2.60689941# - 08 
3.31264027# - 06 
9.24689456# - 05 


Figure 6: Graph of equation of Example 4 in comparing the current method and Wazwaz solution 
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4.4.1 Example 5 

According to Eq.(??), if f{x) = l,g(y) = 4(2e y + e%),A = 0 and B = 0 the equation has been defined as 
follows: 


y"(x) H— y'(x) + 4(2e y + e 2 ) = 0,a:^0 
x 


(42) 


with boundary conditions 


3/(0) = 0, ^/(O) = 0 

which has the following analytical solution: 

y{x) = —2 ln(l +x 2 ). (43) 

We have purposed to RCS method to solve ?? therefore, we construct the residual function as follows: 


Res(a;) = ^£ N y{x) + x) + A(2e^ Nv(x) -t-e^ 1 ’). 

dx z x ax 


(44) 


N + 1 points with equal distance from each other in [0, q\ have been substituted in Res{x ): 

Res(xj) = 0, j = 0,1,2,3,..., N 

A non-linear system of equations has been obtained. Some numerical methods should be used as Newton’s 
method. The Maple version 17 has been used to solve the non-linear system of equations. Solving these 
types of equations system, Maple software has been used the Newton’s method and advanced algorithm. 
The Maple /solve has been used for all of unknown coefficients with an initial zero value. The unknown 
coefficients a, will be obtained and approximation of ^nv{x) is the result. 

?? has shown the comparison of y(x ) obtained by the new method proposed in this paper with N=46, and 
the analytic solution?? .The resulting graph of ?? in comparison to the presented method and the analytic 
solution?? are shown in ??. 

Table 11: Comparison of y{x), between present method and exact solution for Example 5. 


X 

Present method 

Exact value 

Error 

0.00 

0.00000000E + 00 

0.00000000E + 00 

0.00000000E + 00 

0.01 

— 1.999854441? — 04 

— 1.99990000E — 04 

4.55642163E 

-09 

0.10 

— 1.990064341? — 02 

— 1.99006617E — 02 

1.83196007E 

-08 

0.50 

-4.46287089E - 01 

-4.46287103E - 01 

1.39323760E 

-08 

1.00 

— 1.38629436E + 00 

— 1.38629436E + 00 

4.32735003E 

-09 

2.00 

-3.21887583E + 00 

-3.21887582E + 00 

3.97273014E 

-09 

3.00 

—4.60517019E + 00 

-4.60517019E + 00 

5.26229016E 

-09 

4.00 

-5.66642669E + 00 

-5.66642669E + 00 

4.81089035E 

-09 

5.00 

— 6.51619308E + 00 

— 6.51619308E + 00 

3.89665988E 

-09 

6.00 

-7.22183583E + 00 

-7.22183583E + 00 

1.57089008E 

-09 

7.00 

-7.82404602E + 00 

-7.82404601E + 00 

6.89516000E 

-09 

8.00 

-8.34877473E + 00 

-8.34877454E + 00 

1.886484501? 

-07 

9.00 

— 8.81343999E + 00 

-8.81343849E + 00 

1.496180781? 

-06 

10.00 

— 9.23024811E + 00 

-9.23024103E + 00 

7.08015080 E 

-06 
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Figure 7: Graph of equation of Example 5 in comparing the presented method and analytic solution. 

4.4.2 Example 6 

If /( x) = 1, g(y) = — 6 y — 4t/ln(r/), A = 0 and B = 1, one of the Lane-Emden type equations will obtain as 


follows: 


y"(x) + -y'{x) - 6y - 4yln{y) = 0,2 Js 0 


(45) 


with boundary conditions 

2 /( 0 ) = 1 ) 2 /( 0 ) = 0 

which has the following analytical solution: 

y{x) = e x2 . (46) 

We have purposed to RCS method to solve ?? therefore, we construct the residual function as follows: 


Res(a’) = -^- 3 |jv ry{x) + - 6^ N y(x) - 4f w y(a;) ln(| Ar y(a:)). 

ax z x ax 


(47) 


N + 1 points with equal distance from each other in [0, q] have been substituted in Res(x ): 

Res(xj) = 0, j = 0,1,2,3,..., N 

A non-linear system of equations has been obtained. Some numerical methods should be used as Newton’s 
method. The Maple version 17 has been used to solve the non-linear system of equations. Solving these 
types of equations system, Maple software has been used the Newton’s method and advanced algorithm. 
The Maple /solve has been used for all of unknown coefficients with an initial zero value. The unknown 
coefficients at will be obtained and approximation of £a ry(x) is the result. 

?? has shown the comparison of y{x) obtained by the new method proposed in this paper with N=40, and 
the analytic solution?? .The resulting graph of ?? in comparison to the presented method and the analytic 
solution?? are shown in ??. 

The logarithmic graph of the absolute coefficients of Chebysliev of Second kind functions of ?? is shown 
in ??. This graph shows that the new method has a proper convergence rate. 
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Table 12: Comparison of y(x), between present method and exact solution for Example 6. 


X 

Present method 

Exact value 

Error 

0.00 

1 . 00000000 E + 00 

1 . 00000000 £ + 00 

0.0000000015 + 00 

0.01 

1 . 00010001 ^ + 00 

1.0001000LE + 00 

5.0250026415 - 

11 

0.02 

1.00040008^ + 00 

1.00040008.E + 00 

1.5856005415 - 

10 

0.05 

1.00250313^ + 00 

1.0025031315 + 00 

2.5449997715 - 

10 

0.10 

1.01005017E + 00 

1.01005017.E + 00 

2.65950151 E- 

10 

0.20 

1.04081077E + 00 

1.0408107715 + 00 

3.0350011415 - 

10 

0.50 

1.28402542E + 00 

1.2840254215 + 00 

4.3091996 6E - 

10 

0.70 

1.63231622E + 00 

1.6323162215 + 00 

5.87929927 E - 

10 

0.80 

1.89648088^ + 00 

1.8964808815 + 00 

7.86489984 E - 

10 

0.90 

2.24790799E + 00 

2.2479079915 + 00 

9.62370184 E - 

10 

1.00 

2.71828183E + 00 

2.7182818315 + 00 

1.22766020 E - 

09 


Figure 8: 


Graph of equation of Example 6 in comparing the current method and analytic solution. 


Figure 9: Logarithmic graph of absolute coefficients \ai\ of RCS function of Example 6 
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4.4.3 Example 7 

If f(x) = —2(2x 2 + 3),g(y) = y,A = 1 and B = 0, one of the Lane-Emden type equations will obtain as 
follows: 


y"(x) H— y'(x) — 2(2x 2 + 3 )y = 0, x ^ 0 


(48) 


with boundary conditions 

y(0) = l,y'(0)=0 

which has the following analytical solution: 

y{x) = e x \ (49) 

We have purposed to RCS method to solve ?? therefore, we construct the residual function as follows: 


Res(a;) = ^£ N y(x) + --^£ N y{x) - 2(2x 2 + 3 )£ N y{x). 
dx z x ax 


(50) 


N + 1 points with equal distance from each other in [0, q] have been substituted in Res(x ): 

Res(xj) = 0, j = 0,1,2,3,..., N 

A non-linear system of equations has been obtained. Some numerical methods should be used as Newton’s 
method. The Maple version 17 has been used to solve the non-linear system of equations. Solving these 
types of equations system, Maple software has been used the Newton’s method and advanced algorithm. 
The Maple / solve has been used for all of unknown coefficients with an initial zero value. The unknown 
coefficients at will be obtained and approximation of £ nv(x ) is the result. 

?? has shown the comparison of y{x) obtained by the new method proposed in this paper with N=40, and 
the analytic solution?? .The resulting graph of ?? in comparison to the presented method and the analytic 
solution?? are shown in ??. 

The logarithmic graph of the absolute coefficients of Rational Chebyshev of Second kind functions of ?? 
is shown in ??. This graph shows that the new method has a proper convergence rate. 

Table 13: Comparison of y(x), between current method and exact solution for Example 7. 


X 

Present method 

Exact value 

Error 

0.00 

1 . 00000000 E + 00 

1 . 00000000 E + 00 

0 . 000000002 ? + 00 

0.01 

1 . 00010001 E + 00 

1 . 0001000 TE + 00 

4.583000652? - 

11 

0.02 

1.00040008E + 00 

1.00040008.E + 00 

1.072701932? — 

10 

0.05 

1.00250313E + 00 

1.00250313E + 00 

1.383000382? - 

10 

0.10 

1.01005017E + 00 

1.01005017E + 00 

1.27110100 E- 

10 

0.20 

1.04081077E + 00 

1.040810772? + 00 

1.444400162? - 

10 

0.50 

1.28402542E + 00 

1.28402542E + 00 

1.783400092? - 

10 

0.70 

1.63231622E + 00 

1.63231622E + 00 

1.981199652? — 

10 

0.80 

1.89648088E + 00 

1.896480882? + 00 

2.87079915 E - 

10 

0.90 

2.24790799E + 00 

2.24790799 E + 00 

3.109499322? — 

10 

1.00 

2.71828183E + 00 

2.718281832? + 00 

3.63560293 E - 

10 
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Figure 10: Graph of equation of Example 7 in comparing the presented method and analytic solution. 


Figure 11: Logarithmic graph of absolute coefficients \ai\ of RCS function of Example 7 
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5 Conclusions 


The main idea of this paper is to construct an approximation to the solution of nonlinear Lane-Emden 
equations in a semi-infinite interval. A set of Rational Chebyshev of the Second kinds of functions have been 
proposed to provide an effective but simple way to improve the convergence of the solution by the collocation 
method, there is a comparison among the exact solution and the numerical solutions of Horedt 189] and 
the series solutions of Wazwaz [34],Liao[35], Singh et al. [50] and Ramos [41 i and the current work. As a 
result, the present work provides an acceptable approach for Lane-Emden equations. The logarithmic figures 
of absolute coefficients result in the approach which has exponentially convergence rate. Finally, the most 
important concern of spectral methods is the choice of basis functions, there are some properties for basis 
functions such as rapid convergence, easy to compute and completeness, which means that any solution can 
be represented to arbitrarily high accuracy by taking the truncation N are sufficiently large. 
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